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Abstract 

The tumor microenvironment is known to play a key role in altering the properties and behavior of nearby cancer cells. Its 
influence on resistance to endocrine therapy and cancer relapse, however, is poorly understood. Here we investigate the 
interaction of mammary fibroblasts and estrogen receptor-positive breast cancer cells in three-dimensional culture models 
in order to characterize gene expression, cellular changes, and the secreted protein factors involved in the cellular cross-talk. 
We show that fibroblasts, which are the predominant cell type found in the stroma adjacent to the cancer cells in a tumor, 
induce an epithelial-to-mesenchymal transition in the cancer cells, leading to hormone-independent growth, a more 
invasive phenotype, and resistance to endocrine therapy. Here, we applied a label-free chemical imaging modality, Fourier 
transform infrared (FT-IR) spectroscopic imaging, to identify cells that had transitioned to hormone-independent growth. 
Both the molecular and chemical profiles identified here were translated from cell culture to patient samples: a secreted 
protein signature was used to stratify patient populations based on gene expression and FT-IR was used to characterize 
breast tumor patient biopsies. Our findings underscore the role of mammary fibroblasts in promoting aggressiveness and 
endocrine therapy resistance in ER-positive breast cancers and highlight the utility of FT-IR for the further characterization of 
breast cancer samples. 
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Introduction 

More than 70% of breast cancers diagnosed in the US are 
estrogen receptor positive (ER"^) [1,2]. ER"^ tumors generally have 
more favorable prognoses compared to other subtypes and can be 
treated with targeted endocrine therapies such as tamoxifen [3]. 
Though many ER""" patients initially respond favorably to targeted 
therapy, up to 30% of treated cancers recur [3,4]. For patients 
with recurrent disease, the five-year survival rate drops to 20%, 
with a median survival of 12—24 months [5]. Therefore, it would 
be advantageous to identify at the time of initial diagnosis the 
patients who wiU not respond to endocrine therapy in the long- 
term so that their care can be managed differently. The factors 
underlying recurrence arising from endocrine resistance are not 
fully understood, but it is increasingly appreciated that the 
microenvironment of the tumor cells can play a critical role in 
impacting the behavior of the cancer cells [6,7]. 



To understand the molecular factors driving endocrine resis- 
tance and tumor recurrence, we utilized three-dimensional cell co- 
culture models and studied them using molecular profiling and 
chemical imaging. We hypothesized that normal fibroblasts serve 
at the frontline of heterotypic interactions experienced by cancer 
cells because they are the frrst cell type encountered by dysplastic 
epithelium. Further, fibroblasts are encountered in the microen- 
vironment during every stage of disease progression. The 
microenvironment is emerging as a new target for cancer therapies 
[8] . It is now clear that three-dimensional (3D) cultures represent a 
more realistic model for tumors [9,10], and excellent 3D tumor 
models have been proposed [11,12]. However, 3D co-cultures to 
study heterotypic interactions are less widely used [13,14]. Hence, 
we developed and employed a series of 3D co-culture systems to 
investigate the impact of fibroblasts on tumor cell phenotype and 
response to endocrine therapy. Fibroblasts are the most abundant 
cell type in the breast stroma and while they play a role in the 
endocrine regulation of normal breast differentiation, it is not well 
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understood how they affect the response of breast cancer cells to 
targeted endocrine therapy. In order to characterize the influences 
of cancer cell-stromal interactions on therapeutic response, we 
profiled the conditioned medium of the co-culture and defined a 
molecular interaction signature (iSig). The iSig provides mecha- 
nistic insight into tumor progression and the dynamics of cancer 
cell behavior by identifying specific secreted proteins involved in 
cancer cell-stromal cross-talk. When we separated breast cancer 
patient microarray data based on iSig expression levels, we were 
able to predict patient outcome that was comparable to available 
gene expression profiling methods. 

Although uncovering genomic and proteomic dynamics of 
tumor behavior are crucial for understanding the pathophysiology 
of cancer, imaging techniques remain a gold standard of 
determining diagnosis and prognosis in many solid tumors, 
including 1)rcast cancer. Here, we used Fourier Transform 
infrared (FT-IR) spectroscopic imaging [15] for rapid and label- 
free profiling of co-culture samples, integrating molecular infor- 
mation about cellular phenotypes with microscopy. As opposed to 
antibody-based imaging approaches such as immunohistochemis- 
try (IHC) or immunofluorescence, emerging chemical imaging 
technologies use spectroscopy to provide functional and molecular 
information within cells and tissues without staining or the 
requirement for a priori knowledge of molecular transformations. 
By visualizing the chemistry inherent within a sample, cell activity 
changes and disease states in vitro [16,17] and cell types within 
complex tissues [18-20] can be monitored. Although a number of 
studies have reported the development of high performance 
imaging, rapid data processing and classification algorithms, as 
well as applications to histologic analyses of tissues [21-24], few 
studies have related spectral data to underlying molecular 
transformations. 

In this study, we have monitored changes in gene expression 
and cell phenotype associated with cancer cell-fibroblast interac- 
tions and identified chemical changes with FT-IR spectroscopy in 
a "bottom up" approach to identify specific cellular contributions 
of the microenvironment on cancer cell behavior and response to 
endocrine therapy. The first part of this study reports the step- 
wise, validated development of a 3D co-culture model that put 
molecular conclusions in a spatial context. Nc'xt, \vc focus on 
molecular transformations and describe a secreted protein 
interaction signature that effectively predicts clinical outcomes in 
patient tumor samples. Finally, we develop the chemical imaging 
signature of the underlying genomic states in endocrine-responsive 
breast cancer. Together, these highlight a pivotal role for 
mammary fibroblasts located in the tumor cell microenvironment 
in determining the phenotypic properties of ER""" breast cancer 
cells and their sensitivity or resistance to endocrine therapy. 
Further, we show how chemical imaging may be used to rapidly 
identify these cellular phenotypes in heterogeneous samples. 

Results 

3D culture systems were engineered to examine cancer 
cell-mammary fibroblast interactions 

In order to study the interactions between ER"*" breast cancer 
cells and mammary fibroblasts, we utilized four 3D culture models, 
shown schematically in Figure \A and IB. We used a number of 
approaches to translate these findings to human disease states 
(Fig. 1 C and ID). The breast cancer model was comprised of ER"^ 
MCF-7 breast cancer cells grown on Matrigel to form 3D spheroid 
structures, denoted hereon simply as MCF-7. Primary human 
mammary fibroblasts (HMF), isolated from normal breast tissue, 
were commercially obtained. Several co-culture geometries were 



used in order to systematically characterize the mechanism by 
which tumor-adjacent stromal fibroblasts exert influence over 
cancer cell growth. An indirect co-culture, the 'sandwich' (MCF- 
7*), represents the compartmentalization of cancer cells and 
fibroblast-rich stroma, as may be observed in locally confined 
tumors. Because the fibroblasts are grown in a separate coUagen 
layer, the two cell types are able to communicate continuously via 
molecular diffusion of soluble factors while remaining in their 
respective layers [16]. After a prescribed co-culture time, the layers 
are easily separated, and molecular analyses can be performed 
without any cell crossover. The direct, or 'mixed', co-culture 
(MCF-7'^), consists of both MCF-7 and fibroblasts grown together 
on Matrigel, allowing for cell-cell interactions arising from direct 
contact. To understand the transient influence of soluble factors on 
cancer ceU phenotypes, a conditioned medium (CM) culture was 
employed which consists of medium from MCF-7'^ added to the 
MCF-7 ceUs. This medium was subsequently analyzed to identify 
the secreted factors that characterize cancer ceU-fibroblast 
interactions in this model system. These models were used to 
comprehensively interrogate the differential fibroblast-cancer ceU 
interactions that influence cancer ceU behavior and resistance to 
therapy. 

Fibroblasts induce an epithelial-to-mesenchymal 
transition (EMT) in breast cancer cells 

We observed that co-culture with HMF induced an epithelial- 
to-mesenchymal transition (EMT) in MCF-7 breast cancer cells. 
The ceUs grown in the MCF-7''' model displayed increased levels of 
SNAIL (7-fold) and SLUG (25-fold) mRNA as well as a 4-fold 
decrease of E-Cadherin mRNA (Fig. all characteristic of the 
mesenchymal cell phenotype after undergoing EMT [25]. The 
effects observed at 3 days were sustained at longer times, 
indicating that the induced EMT phenotype was not transient. 
This phenotype was accompanied by an increase in mRNA levels 
for growth factors, including amphiregulin (Figure SI in File SI), 
which is often associated with a loss of ER and increased 
inflammatory signaling in breast cancers [26,27]. 

To confirm that the phenotypic changes observed in MCF-7^ 
were due to fibroblast signaling and not dependent on the 
sandwich geometry, we also used the MCF-7'^ culture. Similar to 
MCF-7 and MCF-7^, cells still grew as spheroids in tiie MCF-7'^ 
model, with fibroblasts growing at spheroid perimeters but not 
within them (Figure S2 in File SI). MCF-7^ were grown at two 
seeding compositions, one with a fibroblastxancer ceU ratios of 1:5 
(indicating a 'low density' of fibroblasts) and one with a ratio of 1:3 
(a 'high density' of fibroblasts). In the MCF-7^ cultures, SNAIL 
and SLUG were expressed at much higher levels and E-cadherin 
levels were decreased 2-fold or more compared to MCF-7* 
cultures (Fig. 2B]. A density-dependent effect on expression level of 
these genes was observed between the low and high fibroblast 
ratios. Because the mRNA levels shown represent a mixture of the 
cancer ceUs and fibroblasts, mRNA levels for these markers from 
HMF grown alone are also displayed. Basal levels of SNAIL and 
SLUG in HMF are higher than found in MCF-7 ceUs, but this is 
only approximately 30% the level observed in the co-culture, 
suggesting that regulation of the EMT phenotype is synergistic and 
not dominated by one cell type. Reduction in E-cadherin at the 
mRNA level was also confirmed at the protein level by 
immunohistochemistry (Fig. 2C). 

To demonstrate that the presence of HMF also altered the 
behavior of cancer ceUs, MCF-7 grown in a Transwell co-culture 
with HMF were found to be significantly (p = 0.0045) more 
invasive compared to MCF-7 alone (Fig. 2D). Comparing the 
mixed and sandwich culture models, we observed that increasing 
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Figure 1. Three-dimensional co-culture models and analytical approaches in this study. Schematic of the several three-dimensional co- 
culture models we utilized to study the interactions between iVlCF-7 breast cancer cells and primary mammary fibroblasts. The MCF-7 were grown as 
spheroids in a Matrigel overlay culture. Fibroblasts were incorporated either in a direct-contact mixed culture (MCF-7'"') or in a separate collagen layer 
in a sandwich culture (MCF-7 ). To study the effects of paracrine signaling, conditioned medium (CM) studies were done in which CM was taken from 
the mixed culture and used to treat MCF-7 or normal mammary epithelial cells (HMEC) grown alone. The CM was also profiled using protein arrays to 
obtain the secreted protein interaction signature (iSig). We used gene expression and phenotypic assays to study response to hormone and the 
expression of markers of EMT. This molecular profiling approach was correlated to label-free FT-IR spectroscopic imaging and also gene expression 
from patient samples. 
doi:1 0.1 371 /journal.pone.0096878.g001 



both the co-culture time and the relative number of fibroblasts in 
the co-culture induced similar trends in the expression of EMT 
markers in the cancer cells. These results suggest that fibroblast 
interactions play a significant role in altering the ER"*" cancer cell 
phenotype. 

Co-culture with mammary fibroblasts down-regulates 
ERoc in MCF-7 and drives hormone-independent 
proliferation 

We next examined whether co-culture with fibroblasts might 
alter response of ER"^ cancer cells to hormones. As a prelude to 
using our models for this effort, the response of MCF-7 to estrogen 
was confirmed to be similar in monolayer and in 3D culture 
(Figure S3 in File SI). Interestingly, after co-culture, a sustained 
reduction in the mRNA level of ERot in MCF-7^ was observed 
(Fig. 2iA). Though reduced, the remaining ERa in MCF-7'' cells 
was functional, as shown by induction of the estradiol (E2)- 
responsive genes progesterone receptor (PR) and Ki67 mRNA 
after treatment (Fig. 3-B). The basal levels of Ki67 mRNA and cell 
proliferation were increased in the co-cultures, and after 6 days, 
the response to E2 was greatly diminished, indicating that MCF- 



7^/MCF-7^ proliferation became increasingly hormone-indepen- 
dent (Fig. 35 and 3C). Accompanying this hormone-independent 
growth, the antiestrogen tamoxifen (Tam) was no longer effective 
in inhibiting the growth of MCF-7^ (Fig. 3D). 

Conditioned medium from 3D co-cultures induces EMT 
in breast cancer cells and also in normal mammary 
epithelial cells 

To investigate a possible role for paracrine signaling and soluble 
factors in driving alterations in breast cancer cells, we examined 
gene expression changes in MCF-7 cells treated with conditioned 
medium (CM) (Fig. 4A). Notably, ERot levels were halved in CM- 
treated samples, and EMT markers SNAIL and SLUG were 
upregulated with an almost complete loss of E-Cadherin mRNA 
(Fig. 4^). To examine whether the cancer cell-fibroblast CM might 
also affect normal mammary epithelium, we grew normal primary 
human mammary epithelial cells (HMEC) with CM and 
monitored changes in gene expression. Strikingly, the same 
EMT phenotype was observed although less dramatically than 
in MCF-7 cells (Fig. 4B). Based on these results, we propose that 
fibroblast-cancer cell interactions can induce an EMT in nearby 
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Figure 2. Mammary fibroblasts induce an epithelial-to-mesenchymal transition in ER^ breast cancer cells in 3D culture. (A) Over the 
course of 6 days in the sandwich co-culture, MCF-7^ display increased mRNA levels of EMT markers SNAIL and SLUG with a decrease in E-cadherin 
mRNA. (B) Mixed co-cultures were prepared at two fibroblast seeding densities, low- and high-density (relative to MCF-7) and EMT markers were 
modulated up or down as seen with the sandwich co-culture, and a dose-dependent response to the fibroblast presence was observed. (C) 
Immunohistochemistry was used to confirm decrease in overall E-cadherin protein level InMCF-?'^ (D) Co-culture with human mammary fibroblasts 
(HMF) increases the invasiveness of MCF-7 breast cancer cells. 
doi:1 0.1 371/journal.pone.0096878.g002 



normal breast epithelial as well as cancer cells via soluble protein 
signaling. These observations imply that the environment gener- 
ated by fibroblast-cancer cell interactions may promote aggressive 
cell phenotypes in both normal and cancerous epithelial cells. 

After establishing that fibroblast-cancer cell interactions occur 
via paracrine signaling, we analyzed the CM with a protein array 
in order to elucidate the molecular components involved in this 
interaction. Of the 507 proteins measured on the array, 46 were 
found to be upregulated five-fold or greater in the co-culture CM 
compared to naive MCF-7 CM (Table 1). Notably, many of the 
upregulated proteins were members of the matrix metalloprotei- 
nase family (MMPs), which has been previously implicated in 



aggressiveness of breast cancer cells [28] and human tumors [29]. 
Other identified proteins included growth factors and cytokines 
that play a role in recruiting cells involved in the inflammatory 
response. Functional classification of the secreted proteins also 
showed enrichment in components of chemotaxis and ECM 
remodeling. 

Patients with Piigh expression of the secreted protein 
signature are more lil<ely to have poor clinical outcomes 

To examine whether our in vitro secreted protein signature could 
be used to complement currendy-avaHable approaches for 
stratification of breast cancer patient samples, we utilized a 
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Figure 3. Co-culture with human mammary fibroblasts down-regulates ERa and fosters hormone-independent growth of breast 
cancer cells. (A) Fibroblasts down-regulate ERa mRNA levels at 3 and 6 days and this effect increases with longer co-culture time. (B) HiVIF reduce 
the response of MCF-7 cells to estradiol, as measured by PR, and HMF increases basal Ki67 mRNA at 3 and 6 days and eliminates further increase by 
E2. (C) Co-culture with HMF increases the basal proliferation rate of MCF-7, and decreases the fold change of E2 stimulation of proliferation monitored 
at 3 days. (D) MCF-7'"' proliferation fails to respond to increasing concentrations of E2 (left panel) and MCF-7" growth is not inhibited by treatment 
with tamoxifen (right). Response in MCF-7 cells alone is shown for comparison. 
doi:1 0.1 371 /journal.pone.0096878.g003 



publicly available dataset [30] comprised of 109 patients 
diagnosed with Invasive Breast Cancer (IDC), pure Ductal 
Carcinoma in Situ (DCIS), or a mixed diagnosis (IDC with a 
DCIS component). Genes associated with upregulated CM 
proteins (iSig) were used to cluster microarray data from these 
primary tumors. Three primary clusters were found (Fig. 5^). The 
red cluster, characterized by elevated levels of MMPs and TGFP, 
showed a significantly worse prognosis (Fig. bB). Then, ONCO- 
MINE analysis was used to validate results across additional breast 
cancer datasets [31-34]. The iSig was strongly correlated with 
invasion and poor prognosis, and was associated with stromal gene 
profiles in microdissected samples (Fig. 5Cand Figures S5 and S6 
in File SI). The utility of a signature derived from cell culture 
models to stratify patient data based on gene expression suggests 
that, after further work and validation, bottom-up protein 
signatures may be useful to monitor the dynamic states of tumor 
progression. To our knowledge, this is the first study that utilizes a 
secreted protein signature to stratily breast cancer patient 
outcome. 

FT-IR spectroscopic imaging can be used to monitor 
hormonal response of cancer cells 

Our data show that ER level and activity can be dynamically 
altered and that ER-positive cells may not respond to endocrine 



therapy if they have undergone EMT. Given the spatially- 
heterogeneous nature of tumors, an imaging technique to record 
the individual cellular states within a tumor is needed. While IHC 
can provide information about protein expression, it is time- 
consuming and requires antibodies [35]. We hypothesized that 
using FT-IR microscopy could be used to measure the chemical 
response of breast cancer cells to hormone treatment or to 
endocrine therapy without loss of spatial information. We used the 
same 3D culture samples and conditions so that mRNA expression 
and FT-IR signatures could be directiy compared. After Ej 
treatment, we observed a large increase in the C-H vibrational 
region of the spectrum (Fig. 6^), commonly associated with the 
CH2 bending of fatty acids and lipids [36]. These three peaks 
(2960 cm ',2930 cm and 2850 cm ') are also correlated with 
fatty acyl chain peroxidation*'^. Similarly, the peak at 1080 cm 
associated with the phosphate backbone of nucleic acids and 
altered metabohc activity, is increased upon Eg treatment. Thus, 
changes in absorption pattern can be correlated with biological 
changes occurring in ER^ breast cancer cells in response to 
treatment with Eg. 

Next we sought to differentiate between cells having functionally 
active ER and those that do not. Samples from the mixed co- 
culture were imaged, and the average spectra upon treatment with 
Eg are shown (Fig. 6^). Analysis of the spectra revealed that there 
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Figure 4. Treatment of cancerous and normal breast cells with conditioned medium (CM) alters gene expression. (A) Gene expression 
changes in MCF-7 cells treated with CM reveal a 10-fold increase in SNAIL, a 7-fold increase in SLUG, and elimination of detectable E-Cadherin mRNA. 
The increase in EMI markers is accompanied by a down-regulation of ERa mRNA. (B) After treatment of HMEC grown in monolayer culture with CM, 
there were changes in expression of EMI markers and increases in IL-1 p and the growth factor amphiregulin (AREG), suggesting an important role for 
paracrine signaling between the tumor microenvironment and adjacent normal epithelium. 
doi:1 0.1 371 /journal.pone.0096878.g004 



was a loss of the chemical signature associated with E2 response in 
the co-cultured samples, corresponding to the reduction of ERa 
and ERot-mediated response in MCF-7". The fibroblast dose- 
response seen in mRNA levels was also apparent in IR spectra 
(Figure S7 in File SI). The loss of E2 response was observed in the 
IR spectra of CM-treated samples, but to a lesser extent than in 
the direct co-culture (Figure S8 in File SI). The chemical signature 
of Ej treatment and its loss after co-culture with HMF in 3D 
culture reveals that FT-IR imaging may be useful in determining 
ER function in cells in vitro. 

The cell response to tamoxifen also had a distinct associated 
spectral signature (Figure 6B). There was a change in the peaks 
located in the C-H stretching region of the spectrum in MCF-7 
after treatment with Tarn. Additionally, the anti-proliferative 
effects of Tam were specifically seen in reductions in absorbance at 
the 1080 cm peak. Similarly, the loss of response to Tam in the 
MCF-7'^ co-culture is correlated specifically to a spectral 
signature. While the 1080 cm ' peak was substantially decreased 
in MCF-7 after Tam treatment, the peak was unchanged in MCF- 
7^. Thus, this label-free approach can discriminate between cells 
that are responsive and those that are resistant to tamoxifen 
treatment. 

Next, unstained human breast tumor tissue microarrays (TMAs) 
containing samples from 1 48 patients were imaged with FT-IR in 
order to establish proof of principle on whether spectral differences 
in ER status would be apparent in clinical samples. ER expression 
was determined by IHC staining for ERot on an adjacent TMA 
section. Spectra of epithelium from IDC samples with high 
expression of ER (>80% of cells ER""^ were extracted and 
compared to IDC samples with low expression of ER (<20% of 
cells ER^ (Fig 7). A Bayesiaii-type classifier was developed for the 
histologic discrimination of breast tissue (Fig 7 A) using previously 



reported methods [18]. The classifier was used to extract pixels 
corresponding to epithelium or fibroblast regions from breast 
cancer samples with varying levels of ER. Spectroscopic signatures 
of different cell types, namely epithelial/cancer cells and 
fibroblasts, were compared. Several spectral markers could be 
translated from these 3D culture models to human tissue. In the 
epithelium (Fig IB), spectra from samples with high ER expression 
demonstrated similar absorbance in the C-H stretching region 
(3000 - 2825 cm"') to MCF-7 3D spheroid cultures, and the 
overall absorbance in this region was higher than samples with low 
ER expression. Strikingly, significant spectral differences were 
observed in fibroblasts from tumors with different ER levels 
(Fig 7C). These results suggest that IR might be used to identify 
features in the stromal compartment that are associated with 
endocrine response. Spectroscopic differences between human 
breast tumor samples with varying levels of ER expression suggest 
that FT-IR imaging can complement IHC analysis by providing 
additional chemical information about the function of this nuclear 
hormone receptor in breast tumors. 

Discussion 

Our findings reveal that the interaction of ER"** breast cancer 
cells with mammary fibroblasts promotes a more aggressive and 
endocrine therapy-resistant state. Mammary fibroblasts induce a 
downregulation in the level of cellular ERot and promote a more 
rapidly proliferating hormone-independent state in ER"*" breast 
cancer cells. Importandy, the cross-over to hormone-independent 
growth is associated with acquisition of an EMT phenotype, which 
has been previously correlated to invasive disease and poor 
prognosis [37-40,25]. EMT was concurrent with the disruption of 
E2/ERo( signaling and the occurrence of a more invasive 
phenotype. We also found ERa to be downregulated such that it 
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was no longer the prime determinant of tumor growth, 
underscoring this as a potential mechanism by which tumors 
become resistant to endocrine therapy. Of interest, a recent report 
has shown a similar association between low ERa and EMT and 
poor prognosis in human uterine endometrial carcinoma [40]. 
The therapeutic implication is that although women may be 
diagnosed as ER"*", some wiU not respond favorably to endocrine 
therapies because growth of their tumor is regulated by the nature 
of their tumor microenvironment. Therefore, it may be advanta- 
geous to simultaneously inhibit paracrine factors produced by a 
fibroblast-enriched tumor microenvironment in combination with 
the use of anti-estrogen or other endocrine therapies. 

In analyzing fibroblast-cancer cell interactions for molecular 
targets using protein microarrays, we found this interaction to 
result in the expression of a set of regulatory proteins, including M- 
CSF, SDF-1, and MMP family proteins that have been 
individually associated with cancer activated fibroblasts and 
invasive disease [28,41,42,29]. Our time-course experiments show 
that these represent early signaling events in the fibroblast-cancer 
cell relationship and suggest that it is the combined external 
signaling of both fibroblasts and cancer cells that contributes to the 
extracellular milieu that drives this EMT phenotype. Notably, 
while the presence of fibroblasts increased tumor cell invasiveness, 
the proteins secreted into conditioned medium after co-culture 
also greatly affected the properties of normal mammary epithelial 
cells. These observations imply that fibroblast-cancer cell interac- 
tions may put at risk nearby normal mammary epithelial cells, 
initiating early cancer-like phenotypic changes that may be 
relevant to acceleration of early stages of breast tumorigenesis. 

Table 1. Secreted proteins that were increased 5-fold or 
higher in co-culture conditioned medium (CM) compared to 
CM from naive MCF-7. 
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TGF-p Rll 
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Thymopoietin 


MFG-E8 


TIMP-3 
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TIMP-4 


MMP-1 


TRAIL R1/DR4/TNFRSF10A 


MMP-3 


Trail R3/TNFRSF10C 


MMP-1 3 


VEGF R2 (KDR) 


MMP-14 


VEGF-D 


MMP-1 9 


WIF-1 


MMP-25/MT6-MMP 


XEDAR 
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Our observations are supported by work in mouse mammary 
gland showing the impact of mammary stromal cells in changing 
motility and promoting the initiation of cancer-like properties of 
normal mammary epithelial cells [42] . A recent report shows that 
HGF signaling from normal mammary fibroblasts is able to 
promote basal-like breast cancer phenotypes in a 3D culture 
model [43,44]. Further, fibroblasts may be implicated in the 
etiology of diseases such as benign prostatic hyperplasia that have 
been associated with EMT [45] . While many studies have shown 
that cancer-activated fibroblasts (CAFs) alter tumor behavior [46— 
48], our model system captures early stages of tumor-stroma 
interactions. The hallmark of CAFs is stable expression of Ql-SMA, 
and there is some but very mild staining in the MCF-?''' cultures 
(Figure S9 in File SI), suggesting that the fibroblasts used in this 
study are also being altered by the cancer cells. 

Because the secreted proteins identified in these experiments are 
potentially clinically actionable factors, we used the protein 
expression data to obtain a corresponding panel of genes (iSig). 
This in vitro signature could be used to cluster patient samples and 
predict prognosis. Our approach is complementary to existing 
subtyping based on genetic or histologic means and can be 
employed to broaden and enhance the present system. Whereas 
several signatures of breast cancer invasiveness and poor prognosis 
have been reported [49,50], these are derived from a top-down 
approach utilizing statistical patterns found in the gene expression 
profiles of patient tumors. By taking the results of controlled in vitro 
experiments and comparing the data to patient genome-wide 
expression, our complementary 'bottom up' approach can be used 
to provide mechanistic insight into cell-cell interactions promoting 
cancer progression. In this way, genetic information from 
estabhshed tumors found in patients is related back to specific 
changes seen under controlled experimental settings. This allows 
for the determination of gene and protein features that may be 
seen during early stages of tumor cell migration into the 
surrounding stromal microenvironment. 

Finally, we utilized FT-IR spectroscopic imaging to determine 
hormone sensitivity in biological samples. This label-free and non- 
destructive technique was used to identify chemical signatures of 
disease states. We correlated gene expression and a loss of E2/ERa 
signaling with decreases in peak height in the C-H stretching 
region of the spectrum in 3D culture cell samples, consistent with 
prior observations that EMT is associated with changes in lipid 
profiles in epithelial cells [51]. We also confirmed that MCF-7 cells 
become resistant to tamoxifen using spectral features, primarily in 
the peaks associated with nucleic acids. These results direcdy 
correlate optical profiles to cellular behaviors and genomics. 
Translating the cell culture results to patient samples, the in vitro 
chemical signature was also found in invasive breast cancer 
biopsies with differing levels of ER expression. These results 
indicate that FT-IR imaging can potentially be useful, alongside 
IHC and molecular marker analyses, to determine ERa function- 
ality in patient tumor specimens. Though other imaging 
techniques such as positron emission tomography (PET) have 
been used for functional imaging studies before and after therapy 
[52], development of a similar approach at the microscopic scale is 
complicated by the need to also appreciate the cell type, 
morphology, and spatial phenotypic heterogeneity to understand 
the tumor and its microenvironment. FT-IR imaging can provide 
microscopic evidence rapidly and is applicable at the time of usual 
post-biopsy diagnoses in pathology. 

In conclusion, by combining molecular profiling with chemical 
imaging, we have demonstrated that mammary tissue fibroblasts 
can alter therapeutic response to anticancer agents and play a 
crucial role in controlling whether ER* breast cancer cells are able 
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Figure 5. Hierarchical clustering of breast cancer patient samples based on the secreted protein signature reveals distinct 
subgroups. (A) Hierarchical clustering of primary tumors reveals three distinct classes (B) Kaplan-Meier method and the log-rank test was used to 
compare the mean survival rates across the three identified classes. Patients in the red class (i.e. with higher gene expression of the signature 
secreted proteins) are significantly associated with a poorer prognosis. (C) ONCOMINE analysis shows that the protein signature is also overexpressed 
in IDC compared with DCIS in an independent dataset. The protein signature is also significantly correlated with the stroma of breast tumors 
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to respond well to hormone or become resistant to endocrine 
therapies. Thus, interactions with the tumor microenvironment 
can result in deregulation of estrogen receptor signaling and 
alterations in cell survival signaling molecules that provide tumors 
with alternative proliferative and survival stimuli [53-55]. It was 
recently shown that histologically normal fibroblasts present in 
negative tumor margins can promote invasive phenotypes in 
cancer cells [56], which supports our conclusions that the 
interaction between fibroblasts and cancer cells alters tumor cell 
behavior and therapeutic response. The pro-inflammatory micro- 
environment that is generated as a result of epithelial-fibroblast 
interactions is sufficient to stimulate EMT in breast cancer cells 
and, also strikingly, in normal mammary epithelium. Notably, the 
molecules identified from this interaction are also upregulated in 
breast cancers from patients with invasive disease and predict a 
poor cKnical outcome [28,26,57,27,58,29,59,25]. While the study 
described here was restricted to ER"^ tumors, we hope that this 
model will be extended to other molecular sub-classes of breast 
tumors and eventually may be used to direct individualized 
therapies based on inclusion of the characteristics of a patient's 
tumor microenvironment. Our findings provide new molecular 
insights, coupled with changes in chemical imaging by FT-IR, into 
the microenvironment regulation of tumor behavior via interac- 
tion between cancer cells, fibroblasts, and nearby normal epithelial 



cells. These findings should aid in the development of specific 
inhibitors of molecular species in this crosstalk that promote 
growth and progression to a more aggressive, therapy-resistant 
state. Further, we have shown the utility in combining chemical 
and molecular imaging approaches in order to aid in cancer 
patient management by providing additional prognostic informa- 
tion at the time of the initial biopsy. 

Materials and Methods 

Cell Culture 

MCF-7 cells were grown as previously described [27]. Primary 
normal human mammary fibroblasts (HMF) were purchased from 
ScienCell and cultured according to their specifications (ScienCeU 
Research Laboratories, Carlsbad, California, USA). Primary 
normal human mammary epithelial cells (HMEC) were also 
purchased from ScienCeU and cultured according to their 
specifications. HMEC were used prior to passage 3 to ensure 
phenotype stability. For sandwich (non-contact) co-culture, MCF- 
7 cells were cultured in hormone-depleted medium (TM, 
comprised of phenol red-free basal MEM containing 5% 
charcoal-dextran stripped calf serum). 3D cultures were prepared 
by spreading a thin layer of growth factor-reduced phenol red-free 
Matrigel (BD Biosciences, San Jose, California, USA) in multi-well 
plates and seeding a low density of MCF-7 on top (for a 12-weU 
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Figure 6. Fourier transform-infrared (FT-IR) spectroscopic imaging can be used to monitor hormone response in cells in culture. (A) 

A change in FT-IR spectroscopic imaging is seen in IVlCF-7 cells treated with estradiol (E2), particularly in the C-H stretching region (3000 - 2750 cm"') 
and in the peak associated with nucleic acids (1080 cm"'). However, when cells are co-cultured with human mammary fibroblasts (HMF), the 
response to hormone is lost. (B) Spectral changes are also seen in MCF-7 cells treated with tamoxifen (Tam). While there is a slight induction in peaks 
associated with the C-H stretching region upon treatment with Tam, the peak associated with nucleic acids, and therefore proliferation, is decreased 
in MCF-7 cells. This is correlated with the anti-proliferative effects of tamoxifen on ER* cells. In samples that have been co-cultured with HMF, this 
change is not seen at the 1080 cm"' peak, corresponding to the endocrine-resistant growth that was seen using proliferation assays. 
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plate, 100 |xl of Matrigel was used and 8,000 cells were seeded per 
well). The cells were seeded in 1 ml of TM containing 2.5% 
dissolved Matrigel. Fibroblasts were embedded in 600 jxl of a 
2 mg/ ml type I collagen gel (Type I collagen from rat tail tendon, 
BD Biosciences, San Jose, California, USA) at a density of 
350,000 cells/ml and prepared in 24-well plates as previously 
described [16]. MCF-7 cells were grown in 3D for 5-7 days in TM 
before adding fibroblast layer. Fibroblasts were grown for two days 
in 3D gels before co-culture. Co-culture ranged from 3-6 days 
before analysis. For the 'direct' co-culture, 100 [ll of Matrigel was 
spread in each well of a 12-well plate and 8,000 MCF-7 cells were 
seeded as before. After 5-7 days of growth in TM, HMF were 
seeded at a ratio of 3:1 or 5:1 (2,700 or 1,600 HMF per weU 



respectively). The medium was then switched to a 3:1 or 5:1 ratio 
of TM: complete FM medium (ScienCell Research Laboratories) 
for the duration of the co-culture. MCF-7 cells were grown alone 
as a control, and their medium was the 1:1 mix for these 
experiments. 

Proliferation Assays 

All proliferation assays were performed using the WST-1 
reagent (Roche Applied Science, Indianapolis, Indiana, USA) as 
described [60] . Absorbance was read at 450 nm on a PerkinElmer 
Victor X MultUabel Plate Reader (PerkinElmer, Waltham, 
Massachusetts, USA). Data were normalized to the vehicle control 
based on cell type (for example, MCF-7 treatments were 
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Figure 7. Spectroscopic signatures determined using 3D co-culture models can be translated to human breast tissue samples. (A) 

Tissue microarray (TMA) biopsy cores (1.5 mm core diameter) were IHC stained for ERa and also imaged using FT-IR imaging (N-H/O-H band at 
3300 cm^' visualized here for clarity). Images classified using Bayesian classifier are displayed as well to highlight the ability of FT-IR to discriminate 
between cell types in complex samples. Scale bar represents 250 |im. (B) Differences between epithelial pixels in patient samples with high (>80%) 
and low (<20%) expression of ERoc can be seen in peaks at 1080 cm"' (phosphate) and 1030 cm"' (glycosidic bonds). Interestingly, there are more 
apparent differences in these peaks when pixels from fibroblasts are analyzed. Full spectrum (3800 - 950 cm~'), C-H stretching region (3000 - 
2750 cm"'), and biochemical fingerprint region (1750 - 950 cm ') are shown. 
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normalized to MCF-7 vehicle control, while 3: 1 treatments were 
normalized to the 3:1 vehicle control) in Excel and plotted in 
GraphPad. All samples were prepared in duplicate. 

Invasion assays 

HMF (1x10^ were seeded on the bottom of modified Boyden 
chambers while MCF-7 (3x10*) were seeded on precoated filters 
(8 |im pore size) after membrane rehydration (BD Biosciences). 
Following incubation for 48 h at 37°C, cells were fixed in 10% 
formalin buffer and stained using crystal violet. Non-invaded cells 



on the surface of the filter were removed using a cotton swab. 
Invasion was quantified as described [60] by counting the number 
of cells that invaded the filter compared to the total seeded 
number. Data are expressed as mean ± SD and P<0.05 was 
assigned as statistically significant by using unpaired one-tailed t 
tests (GraphPad Software, San Diego, CA). 

Quantitative Real-time PCR and Primers 

Total RNA was isolated from cells using TRIzol (Invitrogen, 
Life Technologies, Carlsbad, California, USA). cDNA was 
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prepared by reverse transcription with M-MuLV Reverse 
Transcriptase (New England Biolabs Inc, Ipswich, Massachusetts, 
USA). mRNA expression was measured using quantitative real- 
time PGR on a high-throughput ABI Prism 790()HT real-time 
instrument (Applied Biosystems, Life Technologies, Carlsbad, 
Califomia, USA) using SYBR Green PGR Master Mix (Applied 
Biosystems) as described previously [61]. 

Immunohistochemistry 

For both immunohistochemistry and FT-IR imaging analysis, 
samples were fixed in freshly-prepared 4% paraformaldehyde for 
1 h, mounted in Histogel (Thermo Fisher Scientific Inc., 
Waltham, Massachusetts, USA), and then paraffin-embedded. 
For paraffm embedding, samples mounted in Histogel were 
dehydrated with serial ethanc)l washes (50%, 70%, 80%, 95%, 
100%, 100% for 3 h cadi) followed by three 2 h clearing steps in 
xylenes, and finally four 1 h paraffin infiltration steps (ParaPlast 
Plus, Sigma-Aldrich, St. Louis, Missouri, USA). 

For IHG, samples were sectioned at 10 |J.m onto poly-L-lysine 
coated slides (Thermo Fisher Scientific, Waltham, Massachusetts, 
USA). Once dry, slides were deparaffmized (2x3 min xylene, 
1x3 min 1:1 xylene/100% ethanol, 2x3 min 100% ethanol, 
1x3 min 95% ethanol, 1x3 min 70% ethanol, 1 x3 min 50% 
ethanol, cold tap water). Antigen retrieval was performed by 
boiling the samples in the microwave for 15 min in 10 mM 
sodium citrate buffer (pH 6.0). Slides were blocked in Tris- 
buffered saline (TBS) containing 10% FBS overnight at 4°C. 
Primary antibody was added in TBS containing 1 % BSA for one 
hour at room temperature. The remainder of the protocol was 
performed as specified using a HRP/DAB detection kit (EXPOSE 
Mouse and Rabbit Specific HRP/DAB Detection IHC kit. 
Abeam, Cambridge, Massachusetts, USA). E-cadherin antibody 
was used at a dilution of 1:400 (E-cadherin 24E10 Rabbit mAb, 
Cell Signaling Technology, Danvers, Massachusetts, USA), a- 
SMA antibody was used at a dilution of 1:200 (a-SMA clone 1A4, 
Dako Denmark A/S, Glostrup, Denmark). 

Antibody Arrays and Classification of the Protein 
Signature 

Antibody Array was from RayBiotech, Inc. (L-Series 507: 
RayBio Label-based Human Antibody Array 1 -Membrane). Array 
analyzes 507 human proteins simultaneously in two samples to 
examine differential expression. Conditioned medium was ob- 
tained from mixed (direct contact) co-cultures after 3 days and this 
was compared with conditioned medium from MCF-7 cells grown 
in 3D for the same length of time. Samples were prepared 
according to the manufacturer's specifications. Protein array 
membrane was imag("d with an Image Quant LAS 4010 
Luminescent image analyzer (GE Healthcare, Waukesha, Wis- 
consin, USA). GO and pathway classifications of the protein 
signature were conducted using web-based DAVID Bioinformatics 
Resources database (Figure S4 in File SI) [62] and MetaCore 
(Thomson Reuters, New York, NY, USA). 

FT-IR Imaging and Image Classification 

For cell culture samples, paraffin embedded samples were 
sectioned at 5 |J,m onto MirrlR IR-reflective glass slides (Kevley 
Technologies, Chesterland, Ohio, USA). Once dry, samples were 
deparaffinized in hexanes for 24 h and then dried before imaging. 
AH tissue was commercially-obtained (BR961 and BR 1003 Breast 
Tissue Microarrays, US Biomax, Inc., Rockville, Mar\'land, USA). 
Unstained sections were placed on barium fluoride substrates 
(Intemational Crystal Laboratories, Garfield, New Jersey, USA), 



deparaffinized in hexanes for 24 hours, and completely dried 
before imaging. Adjacent sections of the arrays were stained for 
morphology (Hematoxylin and Eosin) and for ERof expression 
(IHC). For FT-IR imaging, a PerkinElmer Spothght 400 
(PerkinElmer, Waltham, Massachusetts, USA) equipped with a 
thermal source and a raster-scanning linear array detector was 
used with a Germanium ATR imaging accessory. An NB medium 
apodization was applied, a 1 cm s~' mirror speed was used for 
acquisition, and zeropadding was not used. Background scans 
were taken at 8 cm ' spectral resolution. 150x150 |J,m images 
were collected at 8 cm spectral resolution with 16 scans per 
pixel and a 1.56x1.56 |im pixel size. Data were atmospheric- and 
ATR-corrected on the Spotlight, and further processing was done 
using ENVI-IDL. A minimum noise fraction (MNF) algorithm was 
applied on all images in ENVI-IDL [63] to reduce noise in the 
data. Two regions were imaged per sample, and samples were 
prepared in duplicate (for a total of 4 data sets per condition). 
Because each image represented a combination of cells, Matrigel, 
and Histogel, a Bayesian-type classifier was developed as 
previously described [18] that labeled spectral contributions from 
cells within the image. Briefly, 10-15 images of 'pure Matrigel', 
'pure Histogel', and samples containing Histogel, Matrigel, and 
cells were imaged. Several thousand pixels were identified and 
manually labeled as coming from one of three classes: Matrigel, 
Histogel, and cells. A three-class Bayesian-type classifier was 
developed using this training set and was subsequendy validated 
on an independent set of samples. This allowed us to extract pixels 
from the image that corresponded to cells rather than extra[:ellular 
matrix or embedding compounds. AH spectra displayed here 
represent the average spectra from the epithelial cells present in 
the sample. Once spectra were classified and extracted, a baseline 
correction was applied and spectra were further normalized to the 
Amide I peak (1656 cm~') for comparison. Spectra were plotted 
using OriginPro software. Similar approaches to FT-IR imaging 
and classification have been described in the literature 
[22,18,64,65,20,66-68]. 

Data Analysis of Microarrays 

A publicly-available microarray dataset was downloaded from 
GEO (GSE26304) and included 31 pure DCIS patients, 36 IDC 
patients, and 42 patients with 'Mixed' histology [30]. Data were 
transformed into GeneSpring GX 7.3 and chips and genes were 
median normalized and median polished. Hierarchical clustering 
of the genes matching the 46-protein secreted signature was 
performed and displayed using Eisen Cluster and TreeView 
software for analysis and visuahzation. Kaplan-Meier method and 
the log-rank test was used to compare the mean survival rates 
across the groups identified by hierarchical clustering.The 46- 
protein signature was uploaded in ONCOMINE (Compendia 
Bioscience, Ann Arbor, MI), using the following; filters: odds ratio 
>2, and P-value <0.001, epithehal versus stroma or invasive 
cancer versus ductal-carcinoma in situ and cancer type (breast 
cancer). The Oncomine database was also used for visualization of 
significanfly associated datasets. 

Supporting Information 

File SI Figures S1-S9. Figure SI. Amphiregulin mRNA 
expression is increased in MCF-7 breast cancer cells in the 
sandwich co-culture (MCF-7^). Figure S2. Light microscopy image 
showing fibroblasts recruited to MCF-7 spheroids after 3 and 5 
days in culture. Figure S3. Comparison of morphology and 
response to E2 in 2D and 3D MCF-7 models. (A) Micrograph of 
MCF-7 cells grown in 2D and 3D culture. (B) Response of MCF-7 



PLOS ONE I www.plosone.org 



11 



May 2014 I Volume 9 | Issue 5 | e96878 



Integration of Molecular Profiling and Chemical Imaging 



cells to E2 stimulation is similar in 2D and 3D cultures as 
monitored by induction of Eg response genes progesterone 
receptor (PR) and Ki67. Figure S4. Table of DAVID pathway 
association based on iSig. Figure S5. Table of ONCOMINE 
dataset patient characteristics. Figure S6. Expression of iSig is 
correlated with breast cancer stromal expression. Figure S7. FT- 
IR signature of Eg response of MCF-7 cells is based on fibroblast 
density. Figure S8. FT-IR signature of Ej response of MCF-7 cells 
is affected by treatment with CM. Figure S9. Fibroblasts recruited 
to the outside of MCF-7 spheroids display mild a-SMA staining. 
(PPTX) 
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